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It is pointed out that models with condensates have nontrivial renormal- 
ization group flow on the tree level. The infinitesimal form of the tree level 
renormalization group equation is obtained and solved numerically for the 
model in the symmetry broken phase. We find an attractive infrared fixed 
point that eliminates the metastable region and reproduces the Maxwell con- 
struction. 



We have two systematic nonperturbative methods to handle multi-particle or quantum 
systems, the saddle point approximation and the renormalization group. Our goal in this 
letter is to combine these two apparently independent approximation methods in order to 
obtain a better understanding of the instabilities and first order phase transitions. 

We start with the path integral, 
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over the configurations We assume the presence of UV and IR cutoffs thus the dimen- 

sion of the domain of integration is large but finite. 

The saddle point approximation to (|1]) coincides with its perturbative expansion when 
the saddle points are trivial, ip = 0. In this case the elementary excitations, the quasipar- 
ticles, are characterized by the eigenfunctions 0„(a;) of the inverse propagator, G~^{x,y) = 
6'^S[(t>]/6(f){x)6(f){y) . The perturbation expansion is applicable for fixed values of the 

cutoffs if the restoring force of the fluctuations to = iplx) = is nonvanishing, i.e. 
the eigenvalues of the inverse propagator, A„, are positive. When the absolute minimum of 
the action is reached at = iplx) 7^ then we follow the strategy of the saddle point 
expansion and the elementary excitations are the fluctuations around iIj{x) 7^ 0. The saddle 
point expansion is applicable as long as A„ > 0. The zero modes, the elementary excitations 
with A„ = 0, correspond to continuous symmetries and are integrated over exactly. If the 
inverse propagator has too many small eigenvalues the saddle point expansion breaks down 
and strong fluctuations develop around ip as we remove one of the cutoffs. This possibility 
brings us to the renormalization group method which is supposed to deal with such problems 

i- 

The basic idea of the renormalization group is the subsequent integration in (|l]). 



me 



-j-SNm^ (2) 



where the field is expanded according to a suitable chosen basis, = Y.n=i 'Pn^n{x), and 
N < 00 plays the role of the UV cutoff. The effective theory for the first N' modes is defined 
through the effective action S^' by the help of the blocking transformation for AA^ = N — N' 
number of modes, 

= f d(f)^,_^, ... [ rf0^e-^^^['^l. (3) 



The elimination step is usually followed by a rescaling in order to restore the cutoff to 
its original value. This rescaling, together with its result, the anomalous dimension, is an 
important device to distinguish the trivial (tree-level) and the dynamical (loop-level) scale 
dependence. 

The fluctuations of the IR modes with n < N can be described within the effective theory 
given by the action 

SN[(t>\=T.9ms3[<P^d,<P,^<P,..], (4) 

j 

where Sj is a suitable chosen complete set of local functionals without keeping the UV 
modes n > N present explicitly. Thus the IR modes decouple from the UV ones since the 
correlations generated by the latters are contained in the numerical value of the effective 
coupling constant gj{N). This trivial observation gives rise a powerful approximation scheme 
by the truncation of (|) when any IR fleld conflguration can be used to reconstruct the 
same overdetermined effective action. In the inflnitesimal form of the renormalization group 
method, 0, where AA^ << N the small parameter of the loop expansion in (|^) is hAN/N, 
thus we expect that the saddle point approximation is apphcable in (|^). The Wegner- 
Houghton equation is obtained by eliminating the plane waves within the shell k — Ak < 
\p\ < k of the momentum space. 
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dk Ak 2 6(f)6(f) h 6(j) ' \ 6(f)6(f) J ' 6(f) 
V + d.. 5S[<I>] 2-ri-d 6S[cl>] 

where the trace and the "■" operation is taken in the subspace of the ehminable modes. 
The second hne generates the rescahng assuming the anomalous dimension rj (input) in 
dimension d. The second order approximation is used for the action as the functional of 
the eliminable field variables which is applicable only if the saddle point amplitude is 0{h). 
Our interest in this work is the case where there is a saddle point during the elimination 
of the modes and we shall see that there is no gurantee in general that the saddle point 
remains 0{h) during the evolution. Thus we need a more reliable evolution equation which 
is applicable without imposing the condition 0{h) on the saddle point amplitude. 

Systems with instabilities display fluctuations with large amplitudes which can be taken 
into account by means of the saddle point approximation. In order to demonstrate the 
importance of the saddle point during the blocking on the renormalized trajectory we shall 
consider the 0^ model in the symmetry broken phase in a box with linear size L and with 
periodic boundary conditions. The bare action of this model is S'a = / + 
[/a(0)] where Ua{(P) = — m^/20^ + g4>'^/4:\ (A is the UV cutoff) and we search for the 
vacuum in the presence of the constraint L~'^ J d'^x < 4>{x) >= The constrained vacuum 
displays spinodal instability or metastability when it is unstable against fluctuations with 
infinitesimal or finite amplitude, respectively. The spinodal instability can be detected by 
local methods, mode by mode inspection. In fact, each fluctuation of the form 

M^) = ^kc'''^'"'' + ^fce-^^^^> = 2pk cos(A;^x^ + ak) (6) 

with infinitesimal amplitude, pjt ~ 0, represents an independent unstable mode of the tree 
level theory when < m? —g^^"^ /2. The Legendre transform of the tree level theory suggests 
that the vacuua with 2m^ / g < $^ < Grn^/g are metastable. The vacuum is stable when 
$^ > &m? /g. It is difficult to identify the possible instabilities of the true vacuum where 
the local, mode by mode analysis is unreliable. This is because the decay of the unstable 
vacuum udergoes large amplitude modifications either in the metastable or the spinodal 
instable region. 

We show now that the renormalization group can be used to deal with the modes of the 
unstable vacuum in a simple one by one manner. The blocking transformation we employ 
consists of the elimination of the modes with momentum k — Ak < \p\ < k, where k is the 
current cutoff and A A; is a momentum parameter smaller than any characteristic momentum 
scale of the system. The effective coupling constants can be defined in the leading order 
of the gradient expansion by Uk{(p) = Hj gjik) / jlcfP . Since any infrared field configuration 
should yield the same effective action we take the simplest choice, a homogeneous field 

The loop expansion for (j^) yields the general result Sn' = J2j^''S^N'. The perturbation 
expansion retains the term 0{W) with j > and the tree level, 0{?f) contribution represents 
a non-perturbative piece what must be considered before the perturbative pieces are taken 
into account. Since the tree level expressions are different from the loop corrections the tree 
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level renormalization, if exists, consists of essentially new and different expressions than the 
loop contributions which have been considered so far |Q. Our constraint generates non- 
trivial saddle points and we consider here the leading order, tree level renormalization only. 
In this approximation the naive scaling is correct {rj = 0, Z = 1), and we ignore the rescaling 
step of the renormalization group method for the sake of simplicity. 

The saddle point of the blocking transformation can be written in the form ipk{x) = 
Y^' ijjpe^P^ where the prime denotes the summation for k — Ak < \p\ < k. The blocked 
potential is given by 



/ Me-s5'=[*+'']l = min [ / d^x (hd^^Y + Uki<l> + ^) 
J J W U \2 



+ 0{h) (7) 



where the fluctuation r] contains only the modes within the shell [A; — A/c, k]. It is worthwhile 
noting that (|^ reduces to the usual local potential approximation of the Wegner-Houghton 
equation if the saddle point vanishes. There might be several minima if) in which case 
one should sum over them. We retain only the single plane wave saddle points, (^). The 
motivation of this rather drastic approximation is the assumption that the saddle point (P) 
seems to be optimal from the point of view of the energy-entropy ballance. To see this 
first note that when Ak —>■ the saddle point is non-vanishing in the momentum space 
on the sphere with radius k only. As far as the saddle point on this sphere is concerned, 
the introduction of other additional plane waves would increase the kinetic energy. The 
phase ak = a-k is a zero mode for each plane waves, corresponding to the translation 
invariance. The more involved saddle points built by several plane waves have the same 
single translational zero mode. Thus the kinetic energy-entropy balance is optimized for the 
plane waves. The resulting tree-level blocking relation is 



Uk-Ak{^) = min 



1 /"i 

fcV^ + - / du f/fc($ + 2pcos(7rii)) 

2 J-i 

1 /"i 

k^pl + - j ^duUk{<^ + 2pkCos{7ru)). (8) 



We followed the tree level evolution of the potential by performing numerically the min- 
imization in (1) at each blocking step. The most interesting aspect of the result is that, 
starting from k = \/m? — Ak, the saddle points ipk are nonvanishing (and consequently 
the potential receives a nontrivial renormalization) for < $^ < ^l^^{k) = 6(m^ — k'^)/g. 
Thus the mode coupling provided by our successive ehmination method extends the spinodal 
instability over the vacuua which are seen as metastable by the tree level Legendre transfor- 
mation. This can happen because once the tree level contributions are found at modes with 

< for certain values of $^ the saddle point contributions renormalize the potential in 
such a manner that spinodal instability occures for other, larger values of 

Other important results of the numerical analysis are in order: 

(i) the saddle point amplitude is a linear function of the field, 2pk = — $ + ^vac{k); 

(ii) we find the potential 

whenever the saddle points are nontrivial, i.e. in the unstable region, see Fig.l, where the 
evolution of the potential with = 0.1 and g = 0.2 is shown; 
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(iii) the blocking converges as Afc — > despite the apparent absence of A/c in (|); 

(iv) we checked that the previous results are universal with respect to the choice of the 
symetry breaking bare potential. 

It is possible to understand the result (iii) if we write (||) for two subsequent steps and 
take the difference: 

1 /"i 

2kpl + 2k'^pkdkPk + 2 7 1 '^^ '^''^'^ '^P'' ^os(7rM)) 



+ J du cos^TTu) dkPk d^Uk + 2pkCos{iTu)) 



0(AA;)2, (10) 



which shows explicitely that the correction to the potential due to one blocking step is 
proportional to AA;. It is also worth remarking that the potential (^ is a solution of the 
equation ([l0|) , as can easely be verified. 

Few remarks are in order at this point, (a) The disappearance of the metastable re- 
gion is in agreement with the finding of ref. ||^. The saddle points of the static prob- 
lem reflect the spontaneous droplet formation dynamics 0: Suppose that the stable vac- 
uum bubble in a sphere of radius R and created in the false vacuum has the energy 
E{R) = AirR'^a — AirR^AE/S where a and AE stand for the surface tension and the 
(free)energy difference between the two vacua. The critical droplet size, Rc = a/AE, be- 
yond which the droplets grow can be identified by the inverse of the highest momentum of 
the condensate, R^ = 6/ g{^l^^{0) — $^). This relation establishes a connection between 
the static and the dynamical properties, (b) The last term in the right hand side of (^ 
was added by hand to achieve a continuous matching of the partition function at the in- 
stability. As a result, (^ reproduces the Maxwell construction for the effective potential 
K//($) = Uk=o{^) i.e. it remains unchanged (at the tree level) in the stable region, 

> $^^^(0), and turns out to be constant and continuous for $^ < <I>^^^(0). (c) The poten- 
tial (H) contains only one coupling constant, rn?{k) = d'^Uk{0) = —k"^. The corresponding 
dimensionless coupling constant, in^{k) = rn?{k)/k'^ = —1, is trivially renormalization group 
invariant. In other words we recover the usual scaling laws in the stable region, while in the 
unstable region the potential is a fixed point given by (||). (d) Similar tree level potential 
has been found in the N-component 0^ model with smooth cutoff where the would be 
unstable and the stable regime join and the naive metastable region is wiped out by the 
Goldstone modes. While the possibility of having a homogeneous ipli^) essential part 
of the argument, our result shows the more general origin of the potential. Furthermore 
note that the keeping of a single plane wave mode as a saddle point is not justified when 
a smooth cutoff is used. The analytical structure of the loop corrections in the vicinity of 
the instability indicates the same potential as well, |T0[. The availability of such a diverse 



derivations and the microscopic dynamics-independent form suggests a more fundamental 
origin of this potential. We believe that underlying reason of (|^) is the Maxwell construction 
and the actual form can easiest be identified by means of the renormalization group method 
where the dynamics of each mode can be dealt with individually. In fact, one can show 
that the differentiability of the renormalization group flow, i.e. the existence of the beta 
functions, and the nonvanishing of the saddle points as Ak requires the form (^. To 
see this consider the finite difference 



Ak 2L^Ak J '^'^^ ^ ' QL^Ak 
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where we expanded (|^) in the saddle point, ipk{x). The convergence of the left hand side 
requires either the smallness of the saddle point, ipk = 0(A/c^/^), or k'^ + d%Uk{^) = 0{Ak) 
and d^Uk{^) = 0{Ak). Since the saddle point is nonvanishing the form follows. Note 
the essential differences between the tree and the loop level renormalization: The former is 
nonalytic in g and lacks the usual logrithms of the latter, (e) The result (i) can be obtained 
by assuming (|]) in the unstable region and requiring continuity of 9$f/fc('^') in ^■ 

The saddle point approximation includes the softest fluctuations, the zero modes, and 
one can obtain nontrivial contributions to the correlation functions at 0{h^). In fact, let 
us insert the product of field variables in the integrand of (^) and follow the successive 
elimination of the modes in the path integral. The resulting tree level expression for the 
2n-point function on the background field $ is 



JV[k]JV[a]Ur=ii^kM) 
JV[k]JV[a] 

E n ^(PPm + PP(2j~i)) n.udf^^ pl (12) 



p j=i 



where ipkip) = J d^^-e'^^^ipkix), Qd stands for the solid angle, we integrate over the zero modes 
characterized by the unit vector k{k) corresponding for each value of the cutoff, k, and the 
phase a{p). The sum in the second equation is over the permutations P of the field variables 
and A;^($) = — Whenevere we have an odd number of fields or $^ > ^lad^) 

one of the pj is such that p| > m? — g^^/6, G^^^ = 0. The integration over the zero modes 
is reminiscent of the integration over the possible rearrangements of the domain walls in 
the mixed phase and restores the translation invariance of the correlations. The two point 
function G^^\pi,p2) = '^-^^jwjPpi^^P'^ +P2) shows that the Fourier transform of ijjp can be 
interpreted as the domain wall structure for a given choice of the zero modes. 

The method put forward in this letter, the use of the saddle point approximation for the 
blocking transformation, can be used to handle some of the systems where large amplitude, 
inhomogeneous fluctuations are present. We studied here an Euclidean system describing the 
equilibrium situation. Whenever the time dependence far from equilibrium has nontrivial 



semiclassical limit this method can be extended to include the real time dependence, [ITT 
Another natural continuation of this work is the inclusion of the loop corrections in addition 
to the tree level pieces. The correlation functions computed in such a manner include the 



fluctuations in a systematic manner which is an improvement compared to ref. [|T2[ where the 



master equation was used to describe the most probable values of the correlation functions. 
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FIGURES 




FIG. 1. The potential Uk{^) for different values of k showing the RG evolution towards the 
Maxwell effective potential. 
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